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Abstract 

The properties of randomly evolving special trees having defined and 
analyzed already in two earlier papers (arXiv:cond-mat/0205650 and arXiv: 



cond-mat/0211092l have been investigated in the case when the continuous 
time parameter converges to infinity. Equations for generating functions of 
the number of nodes and end-nodes in a stationary (i.e. infinitely old) tree 
have been derived. In order to solve exactly these equations we have chosen 
three different distributions for the number of new nodes v produced by 
one dying node. By using appropriate method we have calculated step-by- 
step the probabilities of finding n = 1,2,... nodes as well as end- nodes in 
a stationary random tree. Analyzing the results of numerical calculations 
we have observed that the qualitative properties of stationary random trees 
depend hardly on the character of distribution of v. The conclusion to be 
correct that in the evolution process the formation of a rod-like stationary 
tree is much more probable than the formation of a tree with many branches. 
We have established that the probability of finding n nodes in a stationary 
tree depends sensitively on the average value of v and has a maximum the 
location of which is increasing with n but remains always smaller than unity. 
This is also true for the end- nodes. 

PACS: 02.50.-r, 02.50.Ey, 05.40.-a 



1 Introduction 

In previous two papers [T] , |2] we defined and analyzed random processes with 
continuous time parameter describing the evolution of special trees consisting 
of living and dead nodes connected by lines. It seems to be appropriate to 



repeat briefly the characteristic features of the evolution process. The initial 
state Sq of the tree corresponds to a single living node called root which at 
the end of its life is capable of producing u = 0,1, . . . new living nodes, and 
after that it becomes immediately dead. If z/ > then the new nodes are 
promptly connected to the dead node and each of them independently of the 
others can evolve further like a root. The random evolution of trees with 
continuous time parameter has not been investigated intensively recently. 
The main interest since the late 1990s has been focussed on the study of 
non-equilibrium networks ,3j occurring in common real world. The evolution 
mechanism of trees with living and dead nodes may be useful in some of 
biological systems where the branching processes are dominant. 

In what follows, we will use notations applied in P and j2]. It seems to 
be useful to cite the basic definitions. The probability to find the number 
z/ of hving nodes produced by one dying precursor equal to j was denoted 
by fj where j & Z. ^ For the generating function as well as the expectation 
value and the variance of u we used the following notations: 



E{z''} = q{z), E{u} = qi and {u} = q2 + qi - qf, 

where 



dz^ 



2 = 1 

are factorial moments of u. It was shown in [1^ that the time dependence 
of the random evolution is determined almost completely by the expectation 
value qi. In accordance to this the evolution was called subcritical if gi < 1, 
critical if gi = 1 and supercritical if gi > 1. 

In the further considerations we are going to use four distributions for the 
random variable u. As shown in ^Ij the equations derived for the first and 
the second moments of the number of nodes are independent of the detailed 
structure of the distribution of v provided that the moments q\ and g2 are 
finite. We called distributions of this type for v arbitrary and used the symbol 
z/ G a for its notation . In many cases it seems to be enough to apply the 
truncated distribution of v. If the possible values of the random variable 
V are 0,1 and 2 with probabilities /o,/i and f^, respectively, then in the 
previous paper [2] the distributions of this type were denoted by G t. Many 
times it is expedient to assume distributions to be completely determined by 
one parameter. As known the geometric and Poisson distributions are such 



is the set of non-negative integers. 
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distributions. In paper [7 we used the symbols G g and z/ G p to identify 
these distributions. 

The distribution function of the hfetime r of a hving node will be supposed 
to be exponential, i.e. V{t < t} = 1 — e~*^*. In order to characterize the 
tree evolution two non-negative integer valued random functions fie{t) and 
/Xrf(i) are introduced: fie{t) is the number of living nodes, while is that 
of dead nodes at t > 0. The total number of nodes at t > is denoted by 

Clearly, the nodes can be sorted into groups according to the number of 
outgoing lines. Following the notation in j2] the number of nodes with k > 
outgoing lines at time instant t > is denoted by fi{t, k). A node not having 
outgoing line is called end-node. It is obvious that an end-node is either live or 
dead. Therefore, the number of end-nodes fi{t, 0) can be written as a sum of 
numbers of living and dead end-nodes, i.e. 0) = fii{t, 0)-\-fidit, 0) = fio{t). 
Since all hving nodes are end-nodes fie{t,0) can be replaced by fie{t). The 
total number of dead nodes fid(t) is given by fidit) = J2T=ol^d(t, k). 

In this paper we are dealing with properties of f^{t), ^i{t) and fioit) when 
t ^ oo. We will call the random trees arising from a single root after elapsing 
infinite time stationary. 

In Section 2 the basic properties of probability distributions of the num- 
ber of nodes, living and end-nodes are investigated when t — > oo. Special 
attention is paid in Section 3 to the effect of distribution law of the num- 
ber of outgoing lines. Three different distributions of u are investigated. In 
order to simplify the notation, indices referring to different distributions of 
u are usually omitted in formulas. Finally, the characteristic properties of 
stationary random trees are summarized in Section 4. 



2 General considerations 

Let us introduce the notion of tree size which is nothing else but the total 
number of nodes fi{t) at time moment t > 0. We want to analyze the 
asymptotic behavior of the tree size, i.e. the behavior of the random function 
n{t) when t —>■ oo. We say the limit random variable 

jl = lim fi{t), 
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exists in the sense that the relation: 

hm V{fi{t) = n\So} = V{fi = n} (1) 

is true for all positive integers n, where So denotes the initial state of the 
tree. A randomly evolving tree is called "very old" when t =^ oo, and a very 
old tree, as mentioned already, will be named stationary random tree. 

It is elementary to prove that if the limit probability V{fl = n} = pn 
exists, then the generating function 

oo 

g{z) = J2Pn^'' (2) 

n=l 

is determined by one of the fixed points of the equation 

giz) = z q[giz)]. (3) 

It can be shown that if gi < 1, then the fixed point to be chosen has to satisfy 
the limit relation 

limg{z) = 1, (4) 

while if gi > 1, then it should have the property 

\im g{z) < 1 (5) 

and independently of qi the equation g{0) = must hold. The relation (jH) 
means that the probability to find stationary tree of finite size is evidently 
1, if <?! < 1, but if qi > 1, then 



\im g{z) = < 1, 

zU ^ — ^ 

n=l 

i.e. the probability to find a stationary tree of infinite size is equal to 

oo 



W 

n=l 

The proof of the proposition is simple. Let us assume that 

oo 
n=0 



is a probability generating function, i.e. g(l) = Yl^=o fn = ^ and /i 7^ 1. We 
need the following theorem: 
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Theorem 1 // = gi < 1, then q{z) > z, W < z < 1; while if 

g'(l) = gi > 1, then there is a point < 2;o < 1 such that q{zQ) = Zq, i.e. 
q{z) > z, if < z < zo and q{z) < z, if zq < z < 1. 

Let us introduce the function ip{z) = q{z) — z. Since q{z) is convex, i.e. 
aU derivatives are positive in the interval [0, 1] it is evident that 

dif{z) 

is a nondecreasing function of 2;, V 2; G [0, 1]. If g'(l) < 1, then ip'{z) < 0, 
i.e. ip'{z) is nondecreasing, negative valued function of 2; in < 2; < 1. Since 
ip{l) = it is obvious that ip{z) < 0, if < 2 < 1, i.e. q{z) < z, and this 
is the first statement of Theorem H If q'{l) > 1, then ip'{z) > 0, and since 
93(1) = the inequality ip{z) < has to be true for all z < 1 lying near 1. 
On the other side ip{0) = /o > 1, hence there should exist one "^zq in (0, 1) 
which satisfies the equation (p{zo) = and that implies the second statement 
of Theorem ^ 

It seems to be important to investigate how the living nodes behave in 
very old, i.e. stationary trees. Intuitively one can say that stationary trees 
arising in subcritical evolution do not contain living nodes and they have 
finite average size. At the same time, it seems to be quite obvious that 
stationary trees originating in supercritical evolution could have living nodes 
with non-zero probability, i.e. they are entities of "eternal life". 

In order to give more precise answer the generating function of the random 
variable 



fii = lim fii{t) 



t— >oo 



should be derived. It is easy to show that 



00 



z 



n=0 



is one of the fixed points of the equation 




(6) 



^Existence of more than one zq is excluded because q{z) is convex. 
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By Theorem 1, equation © has one fixed point if qi < 1, namely gf\z) 
1, V z G [0, 1], and from this it follows that 



1, ifn = 0, 



0, if n > 0, 



i.e. the probability that a subcritical very old tree does not have living node, 
is exactly 1. If qi > 1, then besides gi{z) = 1 equation © has another 
fixed point given by gf\z) < 1, V z G [0, 1], hence the probability V{(li > 
0} = 1 — (72^^(0) should be larger than zero. In other words, a supercritical 
stationary tree may evolve infinitely with certain non-zero probability. ^ 

The expectation value and the standard deviation of the total number of 
nodes can be used to characterize the size of a stationary random tree. From 
(0), one obtains 



E{j2} 



dz 



= -. , if Qi < 1- 



In a supercritical evolution, i.e. when gi > 1, the expectation value E{/i} 
does not exist, but the limit relation 

lim Ei/ift) e-('?i-^)^*} = 

t^oo 9l — 1 

can be simply proved. We note here that in the critical case, i.e. when 
qi = 1, the average tree size becomes infinite linearly, i.e. we have the 
relation limj^oo E{/i(i:)/(5t} = 1. 

The standard deviation of the total number of nodes in stationary random 
trees is given by 

^^^^ = TT^^^sT^' if 91 < ^• 
(1 - qiy/^ 

In a supercritical state the standard deviation does not exist, but it can be 
readily shown that 

lim D{/i(t) e-2('?^-i)Q*} = -^L= ( 1 + , if gi > 1. 

V 91 - 1 V 9i - 1 / 

■^More precise formulation would be the following: in countable set T of subcritical 
stationary trees the measure of subset containing trees with living nodes is zero, while in 
that of supercritical trees the measure of subset consisting of trees with living nodes is 
larger than zero. 
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Finally, we would like to deal briefly with some of the properties of end- 
nodes in stationary random trees. By using equation (14) of 0, it can be 
shown that the generating function 



9o[z) 



E{^^«} = y where jlo = lim 0) 

' ^ t — ►no 



n=0 



is nothing else than one of flxed points of the following simple equation: 

q[g,{z)]=g,{z) + {l-z) U (7) 
Substituting 2; = we obtain 



n=l 



i.e. ^o(O) = 0, hence P{/io = 0} = 0. 

In the sequel we will investigate the properties of stationary trees when 
q{z) is known. In this case we can obtain exact expressions for probabilities 
= n}, = n}, P{/io = n}, Vra = 1, 2, . . . from the corresponding 

generating function. 



3 Known distribution of u 

3.1 Truncated arbitrary distribution of u 

The generating function of v is given by 

q{z) = fo + f,z + f2Z^ = l + qi{z~l) + {z - if 



(8) 



with the restriction for qi and q2 determined by the equality /o + /i + /2 = 1- 
(See Fig. 1 in [2].) By using Eq. Q and applying (jHJ we have 



l + qi{g-l) + -q2 {g - 1)' 



The flxed point of this equation is 

1 



L 



l-qiz- V(l - qizy + 252-2(1 - z) 



(9) 
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is also the generating function of /i. It is an elementary task to show that 



0, if ^ = 0, 

g{z} = 1, if z = 1, and gi < 1 

1 - 2(gi - l)/g2, if z = 1, and qi > 1. 



(10) 




number of nodes 



Figure 1: Probabilities to find n = 1, 2, . . . nodes in a "very old" {t = oo) tree when the 
distribution of ly is truncated arbitrary. The dark and light bars correspond to subcritical 
and slightly supercritical tree evolutions, respectively. 



Since 



n=l 



Pn Z 



performing expansion of generating function into power series of z, we 
can determine the probabilities pn, ^ n & Z easily. After elementary 
calculations we obtain 



Pn 



r(n + 1/2) 
r(n + 2) 



where 



= gi - 92 + v2g2 \/ 1 - gi + -g2, 
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and 



2q2 \ll-qi + -q2, 



W .-^ V r(j-l/2) T{n- J + 1/2) 
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Figure 2: Probabilities to find n = 1, 2, . . . nodes in a "very old" (t — oo) tree when the 
distribution of v is truncated arbitrary and /i =0.1. The dark and light bars correspond 
to subcritical and slightly supercritical tree evolutions, respectively. 

The probabilities pn, = 1, 2, . . . can be seen in Fig. ^ These are 
the probabihties to find n = 1,2, .. . nodes in a "very old" tree developed 
according to the distribution of z/ G a. The dark and light bars correspond to 
subcritical (gi = 0.95, q2 = 0.5) and slightly supercritical (gi = 1.05, q2 = 0.5) 
tree evolutions, respectively. In the last case Woo = 0.2 and it is not surprising 
that the probabilities p„ for finite n are larger in sub- than in supercritical 
trees. 

Fig. 121 shows the dependence of p„ on n in sub- and supercritical evolu- 
tions, respectively. If the difference qi — q2 = fi is small enough, let us say 
0.1, then one observes a special phenomenon, namely, an "oscillation" in the 
dependence of Pn versus n, what is clearly seen in Fig. |21 
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Figure 3: Dependence of probabilities Pn, n = 2, 3 on qi at q2 = 0.5 

In order to explain the origin of "oscillation" we calculated the depen- 
dence of probabilities p2 and ps on qi at q2 = 0.5. In Fig. 01 we can see that 
Ps > P2 in the interval < gi < q^i\2). For the sake of simplicity the upper 
limit q['^\2) will be called "critical gi". 

The reason of oscillation is trivial. Clearly, if qi — q2 = fi = 0, then there 
are no random trees having nodes of even number, i.e. p2j = 0, j = 1, 2, . . ., 
therefore, it should be an interval < /i < f'f\2j) in which 

P2j<P2j+l Vj=l,2,.... 

We calculated the critical values of qi corresponding to f^\2j) at q2 = 
0.4(0.1)1.0. The results are shown in Table 1. As expected q[^\2j) is slightly 
decreasing with j when q2 is fixed. 

At this point it is worthwhile to underline one of the most characteristic 
properties of stationary random trees. As known, all moments of the total 
number of nodes are converging to infinity when the evolution is supercrit- 
ical and t — s> oo. However, the probability to find a "very old" (stationary) 
supercritical tree of finite size is always larger than zero and, therefore, the 
probability to find a supercritical tree of infinite size is Woo = 1 — S^i; and 
is always smaller than one. 

In the case of truncated arbitrary distribution of u one can write Eq. (|10|1 
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Table 1: Critical values of qf\2j) at q2 = 0.4(0.1)1.0. 



92 


P2 = P3 


P4 = P5 


P6 = Pi 


0.4 


0.553 


0.514 


0.491 


0.5 


0.674 


0.634 


0.609 


0.6 


0.789 


0.749 


0.723 


0.7 


0.897 


0.858 


0.833 


0.8 


0.100 


0.963 


0.938 


0.9 


1.098 


1.062 


1.039 


1.0 


1.191 


1.157 


1.135 



in the form: 



oo 



n=l 



1 - 2(gi - l)/g2, ifgi>l 



and so, one can formulate the statement as follows: a supercritical "very 
old" tree may be finite with probability 1 — 2(gi — l)/5'2, and consequently, 
infinite with probability Woo = 2(gi — l)/q2- It is elementary to prove that 
2(gi - l)/g2 < 1, if Qi > 1- 

Let us get some insight into the behavior of living nodes in stationary 
random trees. From Eq. © we obtain for all \z\ > that 

W _ / 1, if gi < 1, 



9l,2 



/0//2 = 1 - 2(1 - gi)/g2, ifgi>l, 



and from this we can conclude ^ that 

° 1 /0//2 = 1 - 2(1 - gi)/g2, if gi > 1. ^ ^ 

It seems to be appropriate to underline the notion of the second half of 
Eq. (fT!?|) . It is expressing that the probability of finding living nodes in 
supercritical stationary random tree is 1 — /0//2, i-e. 

P{fL, > 0} = 2 
(I2 

If (71 < 1, then /o > /2 and hence the probability Pq must have the value 1. If gi > 1, 
then /o < /2 and it means that the probability p\p should be equal to /0//2. 
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As to the probabilistic properties of end-nodes in stationary trees we have 
to go back to Eqs. ((Tj) and (jHI)- It is elementary to show that the generating 
function of fto is nothing else than 



9o[z) 



l-/i-2/2-(l-/i) Wl 



4/2/0 
(1-/1)^ 



(14) 



It is seen immediately that 



lim^o(^) = ^o(l) 



1, if gi < 1, 

1 - 2(gi - l)/g2, ifgi>l, 



i.e. the probability that the number of end-nodes in a supercritical random 
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Figure 4: Dependence of probabilities pn^ on n at 52 = 0.5 

tree is infinite and is given by 2(gi — l)/g'2- At this point one has to remem- 
ber that in Eq. (fT^ it was shown that 2(gi — l)/g2, at the same time, is 
equal to the probability of finding infinite number of nodes in a supercritical 
random tree. At the first sight the result of this comparison is surprising, 
but considering that the statement expresses only the equality of measures 
characterizing sets of nodes and end-nodes in an infinite tree, it is far not 
unexpected. By expanding Qq^z) in power series at 2; = one obtains 



,(0)- 1 , r(n-i/2) (4/0/2)"-^ 

Pn - r^h ^^^^^ (l-/i)2"-i' ^^^^ 
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Figure 5: Dependence of probabilities pn^ on qi at n — 2,3 and 52 = 0.5 

In Fig. E] we can see how the probabihty of finding n end-nodes in a 
stationary random tree depends on n when qi = 0.95, q2 = 0.5 and qi = 
1.05, q2 = 0.5, respectively. It is remarkable that the probabilities are 
decreasing very rapidly with increasing n, i.e. there are very few among the 
stationary trees containing a large number of end-nodes. The dependence 
oi p^2^ and pf^ on qi is shown in Fig. El at ^2 = 0.5. One can observe the 
formation of a maximum in both curves. The sites of maxima are lying in 
the interval << gi < 1. 

Table 2: Probabilities p„ and pi"-* for small n 



n 


1 


2 


3 


4 


Pn 


fo 


/o/i 


fo if? + /0/2) 


fo fi Ul + 3/0/2) 


AO) 


/o/(l - /l) 


/oV2/(i - f,r 


2/oV|/(l - fif 


5/oV|/(l - fiY 



In order to compare the probabilities pn and p„ the difference dn = 
Pn — Pn^ has been calculated. For small values of n one can obtain dn < 0. 
Fig. ini shows the dependence of dn on n in sub- and supercritical stationary 
random trees. For the sake of simple comparison of probabilities Pn and pn"* 
for small n values Table 2 has been compiled. 
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Figure 6: Variation of the difference dn = Pn ~ Pn' with n in the case of sub- and 
supercritical stationary random trees 

3.2 Geometric distribution of u 

Let us investigate now the size distribution of stationary random trees when 
u is of geometric distribution, i.e 



1 



l + qi{l-z) 
It can be immediately seen that the equation 



(16) 



l + q,[l-g{z)] 



has to be solved and 

^1,2(2;) 



2gi 



1 + gi ± 7(1 + qi)^ - 4gi 



are the two solutions. Since gi{0) = 1 + 1/qi > 1 the root 

1 
2^ 



92[z) = g{z) 



14 



4gi z 



(17) 



has to be chosen as probabihty generating function for /i. Clearly, 

1, ifgi<i, 



lira g{z) 

2tl 



1/qi, if gi > 1, 



;i8) 



i.e. the probability to find infinite number of nodes in a supercritical station- 
ary random tree is w^o = 1 — 1 /gi . 

By using the power series of g{z) it can be shown that the probability of 
finding n nodes in a stationary tree is nothing else, but 



Pr, 



1 T{n- 1/2) (4gi 



vn-l 



0f r{n + l) [l + qi] 



2n-l 



(19) 
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Figure 7: Probabilities to find n = 1,2,... nodes in a "very old" {t = oo) tree when 
the distribution of v is geometric. The dark and light bars correspond to subcritical and 
slightly supercritical tree evolutions, respectively. 

Fig. [7| shows the dependence of the probability pn on n. The dark bars 
denote probabilities to find n = 1, 2, . . . nodes in a stationary tree produced 
in a subcritical (gi = 0.9) evolution. At the same time, the light bars are 
related to probabilities that n = 1,2,... nodes can be found in a stationary 
tree which was produced by a strongly supercritical (gi = 1.5) evolution. 
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We underline again that E{/i} exists and can be calculated from Eq. (fTTjl 
only if gi < 1. 

The probability distribution of the number of end-nodes in stationary 
random trees has some new qualitative features in comparison with those 
obtained in the case u E t. By using Eq. ^ and taking into account the 
generating function we obtain 

1 1-z 

9o{z) 



1 + qi[l - go{z)] 1 + 
and it is elementary to show that the solution 

^o(^) = TTTT-, V [^ + qi + ql + Qi z - r(gi, 

zqi[i + qi) 

where 

r{qi,z) = \l {1 + qi + ql + qi - 4gi(l + gi)^ z, 

is such that the requirement (70(1) = 1, if gi < 1 is met. Expanding Qq^z) 
into power series around z = Q one finds ^ that 

{gi + (1 + gi + g?) C,[X{q,) + Y{q,)]} , if n = 1, 

= { ^ (20) 

^tg^i^ {a [X{q,T + y(gi)"] - ZMi)} , if ^ > 1, 



""The aim of the following elementary consideration is to give help to understand how 
the probability has been calculated. r[qi,z) can be rewritten into the form: 



r{qi ,z)^qi \/ {z - Zi)(z - Z2), 

where 



and 



91 V?! 

Introducing the notations X(qi) = 1/zi and Y{(i\) ~ l/z2 one obtains 

9o{z) = } r 1 1 + 9i + + 9iz - qi^/z^ V[l - X{qi)z\[l - Y{qi)z\\ , 

and this formula is used to get the power series of go{z). 



16 



where 

1 r(n-l/2) 
2v/? r(n + l) ' 

while 



+ + (l + v/gr)2' Qi + {I + Qi) (1 - v^)' 

and 

n-l 
k=l 

It is relevant to note that 



oo 



n=l 



1, if qi < 1, 
1/gi, if qi > 1. 



In other words, the probability to find infinite number of end-nodes in a 
supercritical stationary tree is nothing else than Woo = 1 — 
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Figure 8: Probabilities to find n = 1,2,... end-nodes in a stationary tree when the 
distribution of v is geometric. The dark and hght bars correspond to sub- and supercritical 
evolutions, respectively. 

Fig. ISl shows the probabilities of finding n = 1, 2, . . . end-nodes in a sta- 
tionary random tree when G g. We have seen in Fig. |3]a similar bar chart 
but by making a comparison between the two bar charts we can conclude 
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that the ratio of probabihty p\ to ^2 is much larger in the case of e g 
than in that of z/ G t. In other words, if v is of geometric distribution then the 
probabihty of formation of rod hke stationary random trees is significantly 
greater than in the case of z/ e t. 

The structure of stationary random trees can be better visualized by 
giving the formulas p„ and pn^ for n = 1, 2, 3, 4 explicitly. Table 3 contains 
these formulas. 



Table 3: Probabilities pn and pn for n = 1, 2, 3, 4 



n 


1 


2 


Pn 


1/(1 + 


gi/(i + gi)' 


AO) 


(1 + + gi + ql) 


ql{l + qi)/{l + qi + qlf 



n 


3 


4 


Pn 


2gi/(l + gi)' 


5gi/(l + gi)' 


Pn 




q\{i + qi)K{qi)/iX + qi + qiy 



where 

hz{qi) = 1 + 3gi + ql and K{qi) = 1 + 7gi + 13^1^ + Iql + ql 
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Figure 9: Dependence of probabilities p„, n = 2, 3 on gi in the case of G g. 
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In the case of geometric distribution of v the curves p2 and ^3 versus q\ are 
similar to those we have seen in Fig. Although the maxima are appearing at 
smaller values of q\ than in the case of z/ G t. The dependence of probabilities 
p„, n = 2, 3 on gi is presented in Fig. El 
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Figure 10: Variation of the difference d„ = p„ — Pn with n in the case of sub- and 
supercritical stationary random trees when G g. 

It seems to be instructive to show the variation of the probability differ- 
ence dn = Pn — Pn' with Ti in the case of sub- and supercritical stationary 
random trees. In Fig. ^Jone can see that c/^ > when n> 1. 



3.3 Poisson distribution 

In the sequel we would like to discuss briefly the case when v has Poisson 
distribution with parameter gi > 0. As known, the generating function of v 
is given by 

q{z) = Y.{z^} = e-(i-")^i 

and E{z/} = T>'^{v] = qi. 

It can be easily shown that the generating function g{z) = E{z'^} of the 

random variable lim^^oo f^{t) = /i should satisfy the following equation: 

g{z) = z q[g{z)] = z e'"' e^^^("). (21) 

Applying Theorem ^ to Eq. (j?H) one can formulate the following state- 
ment: 
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If z 1 1 and < gi < 1, then the equation g = e~'^^ e'^^^ , where g = 
Uracil g{z) = Yl'^=oPn! ^'^^ '^'^^^V ^^^^ '''oot in the interval [0,1], and that is 
g = gi = 1, while if qi > 1, then besides gi there is another root in [0,1], 
namely g = g2 < 1. 

The consequence of this statement means that if gi > 1, then the prob- 
abihty to find infinite number of nodes in a stationary random tree in the 
case of G p is nothing else than 

oo 

Woo = 1 - 5^2 = 1 - y^Pn- 

n=0 




Figure 11: Illustration of the appearance of roots 52 < 1 when qi = 1.1 and 1.2. 

In Fig. ^2we would like to illustrate the appearance of roots smaller than 
1 of the equation R = e~'^^ e'^^^ — g = 0. We can see, that the equation R = 
has only one trivial root in [0, 1], namely the gi = 1, when < gi < 1. The 
black points are referring to roots due to four different values of qi. 

The probability to find infinite number of nodes in a stationary random 
tree has been calculated, and the Woo versus qi curve is plotted in Fig. [T21 
It is interesting to note that Woo > 0.5 when qi = 1.4. 

Now, we want to calculate the probabilities pn of finding n = 1,2,... 
nodes in a stationary random tree. Expanding the expression 



z exp(— gi) exp 



-qi 



n=l 
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Figure 12: Dependence of probability Woo on gi > 1. 

into power series of 2; at z = we can step-by-step compute the probabilities 
Pn, n = 1, 2, . . ., in accordance with Eq. (j2^- We obtain that 

Pn = Cnqi [e-'^^gi]", (22) 

where C„, n = 1, 2, . . . are positive rational numbers. The first seven of 
them are given in Table 4. It seems to be hardly possible to obtain explicit 
an formula for C„, but it is an easy task to compose an algorithm for its 
computation. 

Table 4: Coefficients C„ in p„ for n = 1, 2, 3, 4, 5 



n 


1 


2 


3 


4 


5 


6 


7 




1 


1 


3/2 


8/3 


125/24 


54/5 


16807/720 



It follows from Eq. immediately that the probability Pn vs. qi has a 
maximum at 

1 

gi = 1 - -. 

n 

The appearance of maxima is well seen in Fig. ^Jin the cases of n = 2, 3, 4, 5. 

Finally, we would like to discuss the properties of end-nodes in stationary 
random trees evolved according to Poisson distribution of u. In order to 
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Figure 13: Dependence of probability p„ on qi for n — 2, 3, 4, 5. 

obtain the probabilities pn^ n = 1,2, ... we have to solve the equation 
exp(-gi) exp [gi go{z)] = g^^z) + {l - z) exp(-gi), 

where 

oo 
n=l 

By using a special expansion procedure we can get the probabilities pf\ 
pf\ . . . step- by-step. Introducing the notations 

u{qi) = (e'^^ - qir\ 
h{qi) = e"' + 2q,, 
h{qi) = e^"' + Se^igi + 6ql, 
i^(q-^) = e39i + 22e2^igi + 58e'^'qf + 24^^, 
4(g^) = e^^i + 52e^'^'qi + 328e^'^'ql + AAAe'^'ql + 120g^ 
the first six probabilities are given by the following Eqs.: 

pf^ = u{q,), pf = lq^Hq^)]^ = ]:q!is{qiMq^)f , 

2 D 

= ^q%iqiMqi)Y, = ^q'AiqiMqi)?, 
= ^g?4(gi)Kgi)]^^ 




22 



1 2 3 4 5 6 



1=0.4 

.Q 

5 0.3 
o 

Q.0.2 

























p 

qi=0.9 
qi=1.5 








] 


1 




























1 















1 2 3 4 5 6 



number of end- nodes 

Figure 14: Probabilities to find n = 1, ... 6 end- nodes in a stationary random tree when 
v e p. 
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Figure 15: Dependence of probabilities pn\ n = 2, 3, 4, 5 on qi when v e p. 

By using these formulas we are able to show the dependence of pn^ on n 
and qi, respectively. We see in Fig. HHthat the evolution process prefers the 
one end-node structures, i.e. stationary random trees in which the probability 
of finding n > 1 end-nodes is very small. The curves pn\ n = 2, 3, 4, 5 versus 
qi plotted in Fig. are clearly demonstrating that the sites of maxima of 
probabilities are shifted to higher values of qi as the number of nodes n is 
increasing. 
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4 Conclusions 



We have investigated the main properties of evolution of special trees when 
the continuous time parameter tends to infinity. A tree of infinite age is called 
stationary or simply "very old". We have stated that if the limit relations 

lim V{ii{t) = n|5o} = Pn < 1, V n = 1, 2, . . . 

are true, then the random function ii{t) giving the number of nodes in a tree 
at t > converges in distribution to a random variable fi which counts the 
number of nodes in a stationary tree. 

For the generating function Ej^'^j = g[z) a simple equation has been 
derived, namely 

g{z) = zq[g{z)\, (a) 

where q{z) is nothing else than the generating function of the number of new 
nodes v produced by one dying node. It has been proved that if v has finite 
first and second factorial moments, i.e. if E{z/} = qi and E{i/(i/ — 1)} = q2 
are finite, then 'Yl°^=iPn = 1) if < 1) and 'Yl^=iPn < 1) if > 1- It means, 
if Q'l ^ 1) then the probability Woo of finding infinite number of nodes in a 
stationary random tree is zero, while if gi > 1, then that is larger than zero, 

namely, = 1 - Y^n=iPn < 1- 

Here, one has to note that the expectation value E{/i} exists only if the 
tree evolution is subcritical, i.e. if gi < 1. In the case of qi > 1, i.e. in 
supercritical evolution it has been shown that 

lim E{fi{t) exp[-(gi - l)Qt]} = qj{q, - 1), 

t— >oo 

while if qi = 1, then lim^^oo E{/i(t)/Qt} = 1. 

It has been shown also that the generating function go{z) of the number 
of end-nodes flo in a stationary random tree has to satisfy the equation 

9o{z) = q[go{z)] - (1 - z)fo, (b) 

which can be used to study end-node properties. When q{z) = El^;*"} is 
known then by using appropriate method for comparison of coefficients of 
in both sides of equations (a) and (b) we could calculate step-by-step the 
probabilities Pi,pf'\ P2,P2'\--- to find n — 1,2,... nodes and end-nodes, 
respectively, in a stationary random tree. 
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For the calculations, we have chosen three different distributions of v. 
The generating functions of these distributions are given by the following 
formulas: 



q{z) 



{ l + q^{z-l) + \q^{ 
l/[l+q^{z-l)l 



if z/ G t, 



if z/ e 



if z/ G p. 



Analyzing the results of numerical calculations the first impression is that 
the qualitative properties of stationary random trees depend hardly on the 
character of distribution of v. We have seen that in all cases the probability 
to find n = 1 node (or end-node) in a stationary tree is significantly larger 
then to find n > 1 nodes. One can conclude that in the evolution process the 
formation of a rod-like stationary random tree is much more probable than 
that with many branches. 

We have found special behavior in the dependence of p„ on n only in the 
case of z/ G t when V{v = 1} = /i is smaller than a critical value. The 
appearance of the oscillation of p„ versus n is consequence of the following 
trivial statement: when /i = then p2j = 0, j = 1, 2, . . . . 

It has been demonstrated that the probabilities pn and versus qi show 
a maximum the location of which is increasing with n but remains always 
smaller than 1. This property is best seen in the case of Poisson distribution 
of V. 
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